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Abstract 

We apply the theory of Takahashi and Tachiki in order to explain theoreti- 
cally the dependence of the upper critical magnetic field of a S/N multilayer 
on the temperature. This problem has been already investigated in the lit- 
erature, but with a use of an unphysical scaling parameter for the coherence 
length. We show explicitely that, in order to describe the data, such an 
unphysical parameter is unnecessary if one takes into account the boundary 
resisitivity of the S/N interface. We obtain a very good agreement with the 
experiments for the multilayer systems Nb/Cu and V/Ag, with various layer 
thicknesses. 

I. INTRODUCTION 

In trying to describe the experimental data for different kinds of multilayers, such as 
Nb/Cu or V/Ag, Koperdraad calculated upper critical magnetic fields versus temperature, 
using the Takahashi-Tachiki theory for infinite multilayersB. He used as fitting parameters 
the bulk critical temperature of the S-layer , the ratio between the densities of states of 
the two materials Ns/N^-, and the two corresponding diffusion constants, Ds and Z^jv- 

In calculating the magnetic field anisotropy, which is the ratio between the parallel and 
perpendicular upper critical magnetic fields Hc2,\\/ Hc2,j_, two choices were possible for the 



diffusion constants, which lead to two solutions, called the first and the second solutionH^. 
For the first solution, the fitted parameters are close to what one knows from the measure- 
ments. However, the dimensional crossover, typical for S/N multilayers, appeared to lie 
at a much higher temperature than the measured one. In the second solution the upper 
parallel critical magnetic field exhibits a dimensional crossover at a lower temperature than 
the experimental one. A characteristic of this type of solution is that the superconductivity 
nucleation point for the parallel magnetic field shifts from the S-layer at low temperatures 
to the N-layer at higher temperatures, which seems unphysical for a S/N multilayer whose 

= 0. Another unphysical aspect is that the fitted critical temperature for the S-layer is 
larger than the one known for the bulk (8.9 K). Moreover, instead of an expected concave 
2D-aspect of the curve at lower temperatures, the calculations lead to a convex type of curve. 
In order to fit the experimentally observed dimensional crossover with the theoretical one, 
Koperdraad introduced a scaling parameter a for the magnetic coherence length. However, 
the physical interpretation for this free parameter remains an open question. 

Looking for a physical factor which can replace the role of the unphysical scaling param- 
eter in fitting the data, AartJl suggested to consider finite samples rather than the infinite 
ones on which Koperdraad did his calculations. In finite samples one has to face surface 
effects. Model calculations done on finite samplesi show that the surface nucleation of the 
superconductivity is more pronounced for multilayers with thinner layers, but it almost dis- 
appears as one increases the thickness of the layers. Since the fitting problem mentioned 
above showed up particularly for thick-layer systems, taking into account surface supercon- 
ductivity does not bring any essential improvement to the already existing results. 

In the present paper we consider the influence of an S /N interface resistivity, in order to 
get rid of the unphysical parameter a. This is in line with experimental evidence that the 
interfaces of artificial multilayers for metals with a different crystal structure such as Nb/Cu 
are quite rough0. Indeed, we find that a finite boundary resistivity (Rb) allows for a good 
fit with the experimental data. 

The paper is organized as follows. In Section || we summarize the theory of Takahashi 



and Tachiki, and we introduce the boundary resistivity by adjusting the boundary conditions. 
We also illustrate the role of the boundary resistivity on the proximity effect. Section ^TT| is 
dedicated to the numerical results and conclusion. 



II. THEORY INCLUDING BOUNDARY RESISTIVITY 

First we summarize the Takahashi- Tachiki theory for S /N multilayers. The theory starts 
from the Gor'kov equationO for the pair potential A(r), with a space-dependent coupling 
constant V{r), 

A(r) = V{r)kTY: I rfVQ^(r, r') A(r'), (1) 

in which the summation runs over the Matsubara frequencies. By averaging over the impu- 
rity configurations and considering the dirty limit, it was shown that the integration kernel 
Qu) obeys a Green's function like equation 

[2\uj\ + L(V)]g.(r, r') = 2txN{v)6{v - r'), (2) 

where 

Hp 

L{V) = -hD{v){V-—A{v)f. (3) 

This result appears to be equivalent with a different approach going back to UsadelS'i. The 
material parameters V^(r), A^(r), and D{v) are the BCS coupling constant, the electronic 
density of states at the Fermi energy and the diffusion coefficient, respectively. In practice, 
they are treated as being a constant in each single layer. At the interfaces de Gennes 
boundary conditions are imposedjll which require the continuity of and D(r)(V — 
^A(r))F(r), where the pair amplitude F{jc) is related to the gap function A(r) through 

A(r) = V{v)F{v). (4) 

Takahashi and Tachiki provide a way of solving Eqs. (|I|) and (Q) by developing the 
kernel (^^^(r, r') and the pair function F{y) in terms of a complete set of eigenfunctions of 



the differential operator ^^(V). These eigenfunctions are labeled by the parameter A, and 
the eigenvalues are e^. They are solution of the eigenvalue problem 

L(V)^A = CA^A, (5) 

subject to de Gennes boundary conditions. The requirement of the existence of a solution 
for Eq. (|l]) leads to the equation 

det \Sxx' - 2vrA;T5] — ^ I^aa'I = 0. (6) 

cu ^\^\ + ^A 

For finite multilayers in vacuum, the de Gennes boundary conditions insure that there is 
no current flow through the interface between the multilayer and the vacuum. These bound- 
ary conditions read D(r)(V — ^A{r))\zF{r) = for the pair amplitude, at the interface 
with the vacuum. As usual for these type of layered systems, the growth direction coincides 
with the z direction. When applied to the eigenfunctions \1/a, they become '^'^^g^'^'^'* = 0, 
where we made use of the gauge A(r) = {Hz, 0, 0) when the magnetic field is apphed parallel 
to the layers, and A(r) = (0, Hx, 0) for the perpendicular magnetic field. In the absence of 
a magnetic field, the solution of Eq. (^) giving the largest value for the critical temperature 
is the physical one. In the presence of a field, solving this equation allows us to derive the 
Hc2{T) curves. The temperature at which Hc2 — ^> is T^. 

A further step in applying the theory of Takahashi and Tachiki is to consider the effect 
of the S/N interface resistivity. In our calculations, we make use of more general boundary 
conditions rather than the de Gennes ones. Such boundary conditions were investigated 
by Kupriyanov and Lukichevlll and according to Golubovlll and KhushainovEl they can be 
written as: 

d d 1 F(r+) F(r') 

The boundary resistivity i?^ is a parameter which characterizes the barrier which electrons 
encounter at the interface. A source of this resistance comes from the mismatch of the Fermi 
(or electronic) levels, lattice structure and lattice constant of the two composite metals. As 



a consequence, Rb reduces the migration of the Cooper pairs from the S-layer to the N-layer, 
by that diminishing the proximity effect. 

As we will illustrate in the following, Rb modifies the critical temperature Tc of the mul- 
tilayer and the magnetic field anisotropy, defined as Hc2,\\/ Hc2,±- By consequence, including 
Rb as a parameter, the two solutions used by Koperdraad have to be reconsidered. It will 
turn out that in using the boundary resistivity as a free parameter, only one solution will 
be possible for the fitting, instead of the two solutions. This solution fits the experimental 
data, without using any other free parameter, such as the scaling-parameter a. 

Let us first consider the situation in which there is no magnetic field applied to the system. 
As mentioned already, a finite boundary resistivity reduces the proximity effect. This leads 
to a higher multilayer critical temperature than in the case of perfect transparency of the 
interfaces. As a consequence, the bulk critical temperature used to fit the multilayer 
critical temperature will be smaller than the one used by Koperdraad. This leads us in a 
good direction, since the previously used was higher than the measured value. 

As an illustration of the influence of -R^ on the proximity effect, we calculate the depen- 
dence of the critical temperature of a multilayer on the thickness of the layers for different 
choices for the boundary resistivity. The results for an 11-layers Nb/Cu system are shown 
in FIG. |l[ First, one notices that as the layer thickness decreases, the multilayer critical 
temperature converges smoothly towards 0, whereas in the thick layer limit, it converges to 
the bulk critical temperature T^. Further, the curves show that below a certain thickness 
of the layers, dcr, the superconductivity disappears. Moreover, this critical thickness dcr 
decreases with the increasing of the boundary resistivity, illustrating the fact that due to 
Rb, the density of Cooper pairs is more localized in the S-layers of the multilayer, so that 
the system becomes a better superconductor. 

We consider now the presence of a magnetic field. When a magnetic field is applied to 
the system perpendicularly to the interfaces, due to the in-plane motion of the Cooper pairs, 
the influence of the boundary resistivity is weak. However, for the magnetic field parallel to 
the interfaces, the picture looks different. In this situation, the Cooper pairs move such that 
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they cross the interface, which means that they experience the influence of the boundary 
resistivity much stronger. In the presence of a boundary resistivity, the diffusion of the 
Cooper pairs from the S-layers into the N-layers is diminished. The proximity effect is 
weaker, leading to a higher critical temperature for the same magnitude of the magnetic 
field. Thus we can conclude that the boundary resistivity increases the anisotropy ratio 

In addition it appears that the dimensional crossover temperature is shifted towards 
higher temperatures. This means that the first solution is not favourable, whereas the 
second solution has chances to be ameliorated. 

In the following section we will take as a starting point the second solution and we will 
present the corrections which are performed in view of a fitting with the experimental data. 

III. RESULTS AND CONCLUSIONS 

Considering the second solution, its inconveniences consist in the fact that at low temper- 
ature the Hc2,\\{Tc) curve is convex, instead of the well known concave square-root behaviour 
for the 2D systems. Besides, at high temperatures the nucleation of the superconductivity 
lies in the N- layer, which is unphysical for such S/N systems. Moreover, a too large ratio 
NjyDj^/NsDs is used in fitting, in order to obtain the corresponding anisotropy. 

All these shortcomings are remedied by considering a finite boundary resistivity. In FIG. 
we show results for Nb(17lA)/Cu(376A) multilayer. The solid curves are obtained by 
accounting for a finite Rb- The dashed curves are Koperdraad's result, which could be 
improved by using a scaling parameter, still lacking a physical interpretation. The perpen- 
dicular field curves are not very sensitive to the change of the parameters. We fitted the 
points Hc2(Tc2) = 0, Hc2,\\(T^'^''^), and Hc2,±iT*) on the measured critical field curves, rather 
than the points H^2{Tc2) = 0, Hc2,\\{T*), and Hc2,±{T*), used by Koperdraad. Here T^^° is 
the temperature where the dimensional crossover occurs on the parallel magnetic field curve, 
and T* corresponds to the experimental point at the lowest temperature. In Tabel we 
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show the data used in our fitting (T^ , Ds, D^, and Rb), compared to the data used by 
Koperdraad {T^'^ , Df, and D^). For example, in fitting the Nb(17lA)/Cu(376A) system, 
we used Ds = 2.4 cm^/s, = 78 cm^/s and Rb = 3.17yur2cm, instead of Ds = 0.65 cm^/s 
and Dn = 138 cm^/s, used by Koperdraad et al. The latter set is rather unrealistic, while 
the first set compares nicely with the diffusion constants used by Biagi et a/.0 The resistivity 
has the same order of magnitude as the resistivity of Nb at 77K, which is p^b = SfiQcm, 
and it is an order of magnitude larger than the Cu value of 0.2 fiQcm. Since the interface 
can be considered as a dirty mixture, the value of Rb looks reasonable. The use of a smaller 
and more realistic ratio N^D^/NsDs, can be explained as follows. In the absence of a 
boundary resistivity, Rb = 0, the anisotropy at a certain temperature T* is directly related 
to the ratio Nj^jDj^f /NsDs- However, the anisotropy increases when one considers a finite 
Rb, so that a smaller ratio N]\fDi\f /NsDs is necessary to fit the anisotropy of the upper 
critical fields. Besides, as one can notice in FIG. this choice of the diffusion constants 
is such that the convex behaviour of Koperdraad's Hc2,\\{T) curve is turned into a concave 
one, characteristic for a 2D system. Furthermore, in our solution the nucleation of the su- 
perconductivity takes place in the S-layer, as one expects from physical reasons. Clearly, a 
good agreement between theory and measurements is obtained. 

In the same way we fitted the data for two other Nb/Cu multilayers, as well as for a 
V/Ag system. The results are shown in FIG. | for Nb(172A)/Cu(333A), in FIG. | for 
Nb(168A)/Cu(147A), and in FIG. | for V(240A)/Ag(480A) multilayer. The experimental 
data are taken from the literature@lli. 

In conclusion, by focusing on a fit at the dimensional crossover temperature and allowing 
for a finite boundary resistivity, the theory describes the experimental data nicely. By that 
the merit of the scaling parameter introduced by Koperdraad et al. is reduced considerably, 
the more so as up to now this parameter was not assigned with any physical interpretation. A 
finite boundary resistivity, on the other hand, is in accordance with experimental evidence^. 
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TABLES 



The system 


Tcm 


L>5[cmVs] 




Rb [/xOcm] 


rf[K] 


D§ [cmVs] 




Nb(17lA)/Cu(376A) 


9.20 


2.4 


78 


3.17 


9.89 


0.65 


138 


Nb(172A)/Cu(333A) 


9.20 


1.23 


69 


2.07 


9.88 


0.64 


180 


Nb(168A)/Cu(147A) 


9.50 


1.45 


73 


2.38 


9.61 


0.58 


231 


V(240A)/Ag(480A) 


5.47 


1.1 


54 


3.52 


5.70 


0.67 


73.4 
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CAPTIONS AND FIGURES 



FIG. 0. The critical temperature for an 11-layers Nb/Cu system, as a function of the 
layer thickness, for different values of the boundary resistivity Rb, measured in fiQcm. 

FIG. ^ The upper parallel and perpendicular magnetic fields for the multilayer 
Nb(17lA)/Cu(376A). The dots denote the experimental points^!. 

FIG. ^ The upper parallel and perpendicular magnetic fields experimental and theo- 
retical curves for the multilayer Nb(172A)/Cu(333A). 

FIG. ^ The upper parallel and perpendicular magnetic fields experiment allil and theo- 
retical curves for the multilayer Nb(168A)/Cu(147A). 

FIG. ^ The upper parallel and perpendicular magnetic fields experimentallii and theo- 
retical curves for the multilayer V(240A)/Ag(480A). 
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FIG. 1. The critical temperature Tc for an 11-layers Nb/Cu system, as a function of the layer 
thicknesses for different values of the boundary resistivity i?^, measured in /i$7cm. 
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FIG. 2. The upper parallel and perpendicular magnetic fields for the multilayer 
Nb(17lA)/Cu(376A). The dots denote the experimental points^!. 
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FIG. 3. The upper parallel and perpendicular magnetic fields experimentalEJ and theoretical 
curves for the multilayer Nb(172A)/Cu(333A). 
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FIG. 4. The upper parallel and perpendicular magnetic fields experimentalEl and theoretical 
curves for the multilayer Nb(168A)/Cu(147A). 




FIG. 5. The upper parallel and perpendicular magnetic fields experimental^! and theoretical 
curves for the multilayer V(240A)/Ag(480A). 
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